ON THE DISTRIBUTION OF FREE PATH LENGTHS 
FOR THE PERIODIC LORENTZ GAS III 
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Abstract. For r € (0,1), let Z r = {x e R 2 1 dist(z, Z 2 ) > r/2} 
and r r (x,v) = inf{< > | x + tv S dZ r }. Let $ r (i) be the proba- 
bility that T r (x, v) > t for x and v uniformly distributed in Z r and 
S 1 respectively. We prove in this paper that 

1 f lf \ ft\dr 2 /l 
lim sup - — r / $ r [-] — = — + 



In el L \r J r ir 2 t \ t 2 



1 f 1/4 ft\dr 2 (1 
liminf - — r / $ r - — = — + O 



o+ |lne|7 e \r J r n 2 t \t 2 

as t — > +oo. This result improves upon the bounds on in 
Bourgain-Golse-Wennberg [Commun. Math. Phys. 190 (1998), 
491-508]. We also discuss the applications of this result in the 
context of kinetic theory. 



1. Statement of the problem and main results 
1.1. The periodic Lorentz gas. Let r G (0, |) and define 

(1.1) Z r = {x e R 2 |dist(x,Z 2 ) > r/2}. 

Consider a point particle moving at speed 1 inside Z r and being specu- 
larly reflected each time it meets the boundary of Z r . Such a dynamical 
system is referred to as "a periodic, two-dimensional Lorentz gas" . (In- 
deed, Lorentz used the methods of kinetic theory to describe the motion 
of electrons in a metal as that of a collisionless gas of point particles 
bouncing on the crystalline structure of atoms in the metal |13j ). 

The "free path length" (or "(forward) exit time") starting from x G 
Z r in the direction v G S 1 " is defined as 

(1.2) r r (x, v) = mt{t > | x + tv G dZ r } , (x, v) G Z r x S 1 . 

For each v = (vi,v 2 ) G S 1 such that Viv 2 ^ and the ratio Vi/v 2 G 
R \ Q, one has r r (x,v) < +oo, since any orbit of a linear flow with 
irrational slope on the 2-torus is dense (see for instance in pQ, — section 
51, corollary 1, p. 287 — for this well-known fact). 

Set Y r = Z r /Z 2 ; since r r (x,v) = T r (x + k,v) for each (x,v) G Z T x S 1 
and k G Z 2 , the function r r can be seen as defined on Y r x S 1 . 

l 
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Figure 1. The Lorentz gas: Z r and the punctured torus Y r 

1.2. Invariant measure for the Lorentz gas. Let 

V r = dxdv — meas(F r x S 1 ) and fi r = —dxdv . 

Thus \x T is a Borelian probability measure on7 r x S 1 . 

On the other hand, the evolution of the Lorentz gas is governed by 
the broken Hamiltonian flow 

(x(t),v(t)) = (v(t),0), whenever x(t) G Y r , 

v(t + ) = v(t~) — 2v{t~) ■ n x {t)n x {t) , whenever x{t) G dY r , 

where n x denotes the inward unit normal to dY r at point x. As can be 
easily checked, the measure \i r is invariant under the flow (|1.3j) . 
Let m G L°°(S l ) such that 

(1.4) m > and / / m{9)d^ r {x, 6) = 1 . 

J JYrXS 1 

The object of interest in the present paper is the distribution of r r 
under mdfi r , ie. the function <3>™ : R + — > [0, 1] defined by 

(1.5) $™{t) = mdfx r - Prob ({(x, v) G Y r x S 1 1 r r (x, v) > t}) . 

More precisely, we are interested in the asymptotic behavior of 
as r — > + and for large values of t — large compared to 1/r. 
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As explained in [Z| pp. 221-222, knowing the distribution of the 
free path lengths in the small r limit has several important implica- 
tions. For instance, it leads to the correct asymptotic model for the 
Boltzmann-Grad limit of the periodic Lorentz gas — which is not gov- 
erned by the linear Boltzmann equation (equation (10) of J3]) : see 
below. Another application bears on the asymptotic behavior of the 
Kolmogorov- Sinai entropy — or equivalently, of the Lyapunov expo- 
nent — of the periodic Lorentz gas in the small r limit. We refer to jH] 
for a survey of the most recent results and open questions on this sub- 
ject, and to jHj which addresses one of these open problems by methods 
similar to those developed here. 

1.3. Main result. The reference jH] established for each m G /^(S 1 ) 
as in (jl.4j) the existence of two positive constants C m and C' m such 
that, for each r G (0, |) and each t > 1/r 

(1.6) ^<^ {t) <9k. 

rt rt 
(see Theorems B and C in jH]). In ()1.6|) . the upper bound was proved 
by an argument based on Fourier series, while the lower bound was 
obtained by a geometric construction exhausting all possible channels, 
ie. infinitely long open strips included in Z r . The validity of the lower 
bound in (jl.6j) was extended to arbitrary space dimensions in |12j . 

Numerical simulations in [T2] suggest the following questions: for 
each m G L 00 ^ 1 ) as in (jl.4j) does one have, for each t large enough 
(say, for each t > 2) 

(1.7) $™ (-) -> A m (t) ase-^0? 



and, if so, does one have, for some constant C > 
(1.8) A m (t) ~ - as t -> +oo ? 

We have not been able to fully answer ()1.7|) . but, were (jl.7|) true, 
our main result in this paper ( Theorem 11.11 below) would answer (|1.8|) 
by giving an explicit value for C. It confirms the numerical results 
obtained in [12j (see figures 4 and 5 there). Throughout the paper, the 
notation for the convergence in the sense of Cesaro is as follows: 

1 r* dr 

C-limsup /(e) = / means that limsup-; / fir) — = / 

e-*o+ e^o+ |hie|J e r 

for some e* > 0. (A similar notation is used for the liminf and the lim 
in the sense of Cesaro). 
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Theorem 1.1. Let m e L°° (S 1 ) satisfying fTfl ). Let t* > V2 and let 

AW) = C4imsup$™ f — J , and if) = C4iminf 
T/ien 

(1.9) A?(f) ~ > and Am ( f ) ~ -L as f + 00 • 
More precisely, 

8\\m\\ roo 
~ £* - 3 

A serious shortcoming of the result above is the need for averaging 
in r before letting r — > + . It seems however that it cannot be avoided, 
at least by using the techniques of the present paper. 

This leads to a natural question, that of the choice of the measure * 
to define the Cesaro mean in Theorem 11.11 The reasons for this choice 
are made clear by following the proof, but we take this opportunity of 
giving an idea of this proof by providing some qualitative argument in 
favor of this choice. 

The proof of Theorem 11.11 is based on comparing the size r of the 
obstacle with the sequence of errors d n in the approximation by con- 
tinued fractions of v%jv\, ie. of the slope of the linear flow (say, in the 
case where < v% < Vi) . It is natural in this context to renormalize the 
problem by applying to a = v 2 /vi the Gauss map T : a: i— >- — [-] - 
we refer to the appendix for more details on these notions. Lemma ITD 
below and especially the formula d n (a) = ad n _i(T(a)) show that the 
exit time problem with slope a and obstacle of size r is mapped to the 
analogous problem with slope Ta and obstacle of size ar. Hence it is 
natural to define the Cesaro average in Theorem 11.11 with the measure 
— which is the scale invariant (Haar) measure of the multiplicative 
group R^. 

Remark 1.1. In fact, one can prove that A™(t*) = A™(£*) for each 
t* large enough (say, for t* > 10). (In other words C-4im Tw0 + (V) 
exists). However, this result requires some significant improvements of 
the method used in the present paper; they will be described in jS]. 

2. A PARTITION OF T 2 

In 1989, R. Thorn posed the following problem: 

• What is the longest orbit of a linear flow with irrational slope 
on a flat torus with a disk removed? 




;i.io) 



A™(t* 



ttH* 
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Figure 2. The punctured torus Y r and the slit S r (v) 

The answer to this question was found by Blank and Krikorian in [2] 
and is summarized as follows. Without loss of generality, assume that 
the linear flow is x \— > x + tv where v = (cos8, sin9) with 9 G (0, j). 
The removed disk of diameter r is replaced by a vertical slit S r (v) of 
length rj cos 9 with the same center (see figure 2). 

Proposition 2.1 (Blank-Krikorian p. 722). Let r G (0, \), 9 G 

(0, f) and f = (cos 0, sin 6 1 ). Assume that tan# zs irrational. There 
exist three positive numbers Ia(t,v), /s(r, v) and lc(r,v) satisfying 

l A (r,v) < l B (r,v) and lc{r,v) = l A (r,v) + l B (r,v) 

and such that, for any orbit 7 of the linear flow x 1— > x+tv in T 2 \S r (v), 
length^) takes one of the three values l A (r,v), Ib{t,v) orlc(r,v). 

Conversely, given any I G {£a(?", v), Zs(r, v ), lc(r, v)}, there exists an 
orbit 7 of the flow x 1— > x + tv of length I. 

Let v be fixed; orbits of length l A (r,v) (resp. of length l B (r,v), 
lc(r,v)) are referred to as orbits of type A (resp. of type B, C) 1 . 

1 On p. 722 of PI, the sentence "If the moving slit first meets the fixed slit at the 
bottom, the A and £?-orbits are reversed" might be the source of a slight ambiguity 
in the definition of the partition above. In the present paper, the orbits of type A 
are the shortest, consistently with the table on p. 726 of [2]. Hence the roles of 
orbits of type A and B cannot be reversed in the present discussion. 
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Figure 3. The partition (Y A (r,v),Y B (r,v),Y c (r,v)) of 
T 2 \ S r (v) 

Proposition 12.11 defines a partition (Ya(t, v ), Ys(r, v), Yc(r, v)) of T 2 \ 
S r (v), where 

Ya(j, v) = {x G Y r \ S r (v) | x belongs to an orbit of type A} 

and Yb(t, v) and Yc(r,v) are similarly defined. Define further 

Sa(^,v) = {y G S r (v) | the f-orbit starting from y is of type A} 

with analogous definitions for Sb(t, v) and Sc(r,v). Clearly Y^(r, v) 
(resp. Yb(t,v), Yc(r,v)) is metrically equivalent to a strip (parallel- 
logram) of length ^(r, v) and width \Sji(r, v)\ (resp. of length l B (r,v), 
lc(r,v) and width \S B (r,v)\, \S B (r,v)\): see figure 3 below. 

Define now A r to be the exit time in the torus with the slit (instead 
of the disk) removed, ie. 

(2.1) X r (z, v) = m£{t > 0\z + tv e S r (v)} for each z G T 2 \ S r {v) . 
For v = (cos 8, sin 9) with 9 G (0, f) as above, define 

(2.2) ^ r (t, v ) = dx- Ymh({z G T 2 \ ^(v) | X r (z, v) > t}) . 

With the partition of T 2 \S' r (f ) in (Ya(t, v), Y B (r, v), Yc(r, v)) , which is 
metrically equivalent to the disjoint union of three strips as represented 
on figure 3, computing ip r (t,v) in terms of the quantities \S A (r, v)\, 
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Figure 4. Graph of 1 1-> ip r (t,v) 

l A (r,v), \S B (r, v)\, l B (r, v), \Sc(r,v)\ and /c( r , becomes an easy task. 
One finds that 

• if < t < l A (r, v), then 

(2.3) i{> r (t,v) = l- tr; 

• if l A (r,v) < t' < Ib{t, v), then 

ip r (t',v) = l-l A (r,v)r 
~ 1 -(t'-l A (r,v))(\S B (r,v)\ + \S c (r,v)\)cos9; 

• if Ib{t,v) < t" < lc(r,v), then 

(2.5) Mt",v) = (l c (r,v) -t")\S c (r,v)\ cosG ; 

• if t > lc(r, v), then 

(2.6) ip r {t,v) = 0. 

The graph of 1 1 — > ip r (t, v) is presented on figure 4 below. 

So far, the distribution ip r (t,v) has been computed in terms of the 
quantities \S A (r,v)\, l A (r,v), \S B {r,v)\, l B {r,v), \S c (r,v)\ and lc{r,v) 
that characterize the partition (Y A (r, v), Y B (r, v), Ya(r, v)) of T 2 \S r (v). 
These quantities can be expressed (see j2j) in terms of the continued 
fraction expansion of tan 9, in the following manner. 

In the following discussion, we freely use the notations recalled in 
the Appendix below. Let a = tan# e (0, 1); the 6*'s for which a E Q 
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form a set of measure and are discarded in the argument below. With 
the sequence of errors d n in the continued fraction expansion of a, we 
consider the following partition of the interval (0, 1): 

(2.7) (0, 1) = |J I n , with I n = [d n , d n _i) . 

Each interval I n is further partitioned into 
(2.8) 

I n = [J In,k , with I nM = [sup(d n , d n _i - kd n ), d n _x - (fc - l)d n ) , 

l<k<a n 

so that eventually we arrive at the following nested partition of (0, 1): 
(2-9) (0,1) = |J |J I n ,k. 

n>l l<k<a n 

Proposition 2.2 (Blank-Krikorian p. 726). Let 9 e (0, f ) and 

r e (0,±). Let v = (cos 9, sin 9) and set a = tan# (it is assumed that 
a ^ QJ. For i? = - j —e, the integers n > 1 and snc/i that 1 < k < a n 
(a n being the n-th term in the continued fraction expansion of a) are 
defined by the fact that R e I n (since the intervals I n ^ form a partition 
of (0,1): see / TOP ). T/jen 

• l A (r,v) = q n and \S A (r,v)\ = R-d n , 

• l B (r,v) = g n _i + kq n and \S B (r,v)\ = R - (cf n _i - fcd n ) ; 

• lc(r,v) = q n -i + (k + l)q n and \S c {r,v)\ = d n - 1 -(k-l)d n - R. 

Using these values, we arrive at the following expression for ip r (t, v), 
whenever R = 6 I nk (see figure 5 below): 

• if < tcos9 < q n , then 

(2.10) ip r {t,v) = 1 - tr , 

• if Qn < tcos9 < q n -i + kq n , then 

(2.11) A(t,v) = 1 - i?g„ - d n (t cos 9 - q n ) , 

• if g n _i + A;g n < tcos9 < q n _i + (A; + l)q n , then 

i2 ^>r(*, «) = l-ityn - dn[?n-l + ~ l)<?n] 

- [i? - (d n _! - (A; - l)cf„)] (t cos 6-q n -x- kq n ) , 

• if tcos9 > q n -\ + (k + l)q n , then 
(2.13) ip r {t,v) = 0. 

Finally, the discussion above leads to the statistics of the exit time 
X r (z,v) defined in ()2.1j) corresponding to the torus T 2 with the slit 
S r (v) removed. In the case of the torus with the disk removed, the 
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Figure 5. Graph of 1 i-> ip r (t, v) for R e I n ,k 

corresponding exit time r r (x,v) defined in f)1.2|) is related to X r (z,v) 
by the obvious inequalities 

(2.14) X r (x, v) — ~ < T r (x,v) < X r (x,v) + | for each x £ Y r \ S r (v) . 
Define 

(2.15) <p r (t,v) = dx- Prob({x eY r \ S r (v) |r r (x,i;) > *}) \ 
because of (|2.14|) . one has 

(2.16) Vr(*+5,v)<0r(t,v)<^r(t-i,v), *>§• 

The remaining part of the paper uses the evaluation of <p r based on this 
inequality together with the formulas (|2~TUJ) . (I2~TTJ) . (17121 and 

for ?/v. 

3. An ergodic theorem 
Given a £ (0, 1) and e £ (0, 1), we define 

(3.1) N(a,e) = inf{n £ N | d n+1 (a) < e} . 

In terms of the partition ()2.7|) of (0,1), N(a,e) can be equivalently 
defined by the condition 

(3.2) e £ ijv( a ,e)+l • 

We start by recalling the following more or less classical lemma. 
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Lemma 3.1. For ae. a G (0, 1), one has 

N{a,e) ~ -i^lne as e -> 0+ . 

Proof. The definition of N(a,e) and the third formula in Lemma 17.11 
imply that 

iV(a,e)-l iV(a,e) 

(3.3) ~^T j a< -lne< -lnT^a, 

i=o j=o 

where T is the Gauss map (|7.1|) . First we prove that N(a,e) — > +oo 
as e — > + ae. in a G (0, 1). Indeed, let C > and let 

£ c = {q6 (0, 1) | N(a, e)<C, for all e G (0, ±)} . 

For all a G i?c and all e G (0, |) 

c 

— In e < — In T^a ; 

3=0 

therefore, for all e G (0, |), 

— \ne-dg — meas (-EcO < (C + 1) / (— In a)rf^(a) < +oo , 

Jo 

since the measure <ig in ()7.2j) is invariant under T. This implies that 

(3.4) for each C > , cfc/ — meas (i?c) = . 
Because N(a, e) is a nonincreasing function of e, 



for each cc G ( £7 m J 

\m>l / 



(3.5) for each a G j [ J E m j . ,V(n. r) — x 

as e — > + and 
(3.6) 



dg - meas | [J £ m J 

\m>l / 



Secondly, by Birkhoff's ergodic theorem, there exists a (ig-negligeable 
set E' such that 

for each a G (-E 1 '^ , — — lnT^a — >■ / (— In a)dg(a) 

N o ^ 

as AT — > +oo, since the Gauss transformation T of (0, 1) is ergodic with 
respect to the invariant measure dg(a) - - see the Appendix below. 
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for each a G ( E' U (J £ m ] 



1 iV(a,e)-l J 

— ; > — In T^a^ l {— \na)dq(a) . 

N(a,e) ^ 

By ()3.3|) . ()3.5|) and (|3.6|) . one finally obtains that, as e — > + , 

- ln ^ . f\ ,__ AJ _,_ A Cg) 

2 In 2 



for each a G ( -E' U I J £ m ] , — — r- — ► [ (—\na)dg(a) 



(replacing the In under the integral sign by its Taylor series at a — 
1). ' ^ □ 

The main result of this section is the following application of the 
Birkhoff ergodic theorem. We shall need the notations below: 

Ao(a, x) = -x- In d N{a>e - x)+1 (a) , 
Ai(a, x) — — x — hidN^ aie -x\(a) . 

Proposition 3.1. Let f be a bounded nonnegative measurable function 
on R 2 . For each x* G R and ae. in a & (0, 1), one has 

1 /" |lne| x/A ,_ A _, NJ _ 12 [ l F(6)d0 



/(Ao(a, x), Ai(a, x))dx 



as e —> 0, where 

r\Ht)\ 

F{0)= / /(| ln(0)| -y,-y)dy. 



Proof. First we decompose the integral 

»|lne| 

/(Ao(a,o;),Ai(a,a:))da; 

o 



a;* 



/(A (a,x), Ai(a,x))dx 

(3.8) f\lne\ 



+ / f(A (a,x),Ai(a,x))dx 

J\ln djv( a , e )(a)| 

,|lnd i+1 (a)| 

+ / /(Ao(a,a:),Ai(a,ar))da; 
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Now observe that, whenever x belongs to the domain of integration of 

-I lnd i+1 (a)| 

f(A (a,x),A 1 (a,x))dx, 

'|lndj(a)| 

in other words, whenever 

| In di (a) | < x < | In c^ +1 (a) | , then N(a, e~ x ) = I . 
This implies that 

'| lnd i+ i(a)| 

f(A (a,x),A 1 (a,x))dx 

'I In <*!(<*)! 

- In d l+1 (a) 



(3.9) 



f(—x — ln<i; + i(a), —x — \ndi(a))dx 

\ndi(a) 

In^oOM+iM) 

f(\n(di(a)/di + i(a)) - y, -y)dy 



= F(T l 6) . 

Likewise, since / is nonnegative 
(3.10) 

Hlne| 

0< / /(A (a, x),Ai(a, x))dx 

J\ lnd N( a , e )(°<)\ 

/•|lndjv(a,e)+l(°!)l 

< / f(A (a,x),A 1 (a,x))dx 

J I Induce) (") I 

< i^T^^a) . 

Using (|3IBD and flUU) leads to 
(3.11) 

.|lne| 

/(A (a, x),Ai(a, x))dx 



In el 



N(a,e) 
| rnej I N(a 



N(a,e) \ 



1 '° 



In el 



/(A (o!,x), Ax(a,x))dx 



1 



[lH6| 



+ - — -l f(A (a,x),A 1 (a,x))dx. 

I ine l ■/|lndAr(a,e)(a)l 

Observe first that, by its definition 

< F(6) < ||/|| L oo|ln^| which implies that F E L 1 ((0,l),dg) . 
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Birkhoff 's ergodic theorem applied to the Gauss transformation T to- 
gether with Lemma 13. II shows that 

JV(a,e) 



£ F{T l 9) - — / F{9)dg{9) 

l—n J 



N(a,e, 
K ' ' 1=0 

for ae. a G (0, 1) as e — > 0, so that the first term in the right hand side 
of (j3.11|) converges to 

i| f 1 F(6)dg(e) 
Jo 

as e — > + ae. in a. The second term in the right hand side of (j3.11|) 
obviously vanishes ae. in a G (0, 1) as e — > 0. As for the third term, 
because of (13.101) . one has 



I r\ ln t\ 

0<- — r/ /(A (a,x),Ai(a,a;))da; 



~ 'lne| iV(a,e) 



and the right hand side of the above inequality converges to ae. in 
a G (0, 1) as e — > + by Birkoff's ergodic theorem and Lemma 13. 11 □ 

4. Proof of theorem 11.11 

4.1. Step 1: pointwise estimates. The discussion in section El made 
it clear that the natural objects for applying Birkhoff 's ergodic theorem 
are functions involving a finite number (two in this case) of the d n 's 
with n — > +oo. In this first step, we shall reduce the function ip r 
to an expression of this form modulo terms that are small in some 
appropriate sense in the asymptotic regime that we consider — ie. as 
t* -> +oo. 

Throughout this subsection, we set r G (0, |) and 6 G (0, f ) such 
that a = tan# G' Q; let then v = (cos 9, sin #) and i? = —7^. We also 
consider n > 1 and G {1, . . . , a n } (where a = [ai, a 2 . . .]) such that 
i? G J n) fc — this defines n and A; in a unique way, since the I n .kS form 
a partition of (0, 1). 

Lemma 4.1. Under these conditions on a, n k and R 
• the integer k is given by the formula 

(4.1) k = M{1 G N* I rf n _! - ld n <R}= dn ~ l ~ R 

CLri 
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• the denominator q n of the n-th convergent of a satisfies the 
estimate 

1 1 
(4 ' 2) R+(k + l)d n <qn< R+{k-l)d n ' 

• finally, one has 

1 2 
(4.3) q n d n < - , < 1 - q n dn-\ < 



k ' ™ n ~ L k + 1 

Proof. Whenever d n < R < d n _i - - ie. whenever R G J n - - the 
condition i? G amounts to defining k by the formula 



k = inf{Z G N* | rf n _i - ld n < R} 



d n _i — i? 

(If, 



+ 1. 



Next we estimate g n . Since -R G J ni fc, one has in particular the 
inequalities 

(4.4) d n _i -kd n <R< rf n _i - d n (k - 1) , 

which imply 

1 1 

> 1 > <?n , 



+ — 1) d n -i 

by the second inequality in (|7.9j) . This is exactly the upper bound in 
(O- By (JZZZj), 

1 - Rq n - d n (k - l)q n = d n q n „i + d n -tq n - Rq n - d n (k - l)q n 

< q n {d n + dn-i -R - (k- l)d n ) 

< q n {d n + d n -i - d n -i + kd n - (k - l)d n ) 
= 1d n q n 

which gives the lower bound in (|4.2j) . 

The upper bound in (|4.2j) and the fact that R> d n (since R G 
imply the first inequality in ([4.3)1 . As for the second inequality there, 
observe that q n d n -i < 1 by (|7.9|) which establishes the lower bound, 
while the upper bound follows from ()4.2|) in the following manner: 

1 R + d n (k + 1) - d n _i 

q n < 



d n -i d n -± (R + d„(k+l)) 

< . ,„ 2d : < 



d n _ x (R + d n (k + 1)) d n ^(k + l) 
(where the penultimate inequality above follows from (|4.4)l ). □ 
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Lemma 4.2. Let v, a, n, k and R be chosen as above, and let t* > 2. 
Define 



(4.5) 



t* 



l 



R 

d n - 



then for ^ > q n we have 



(4.6) 



ijj r | —,r 



Xr 



t* 



-,v 



R 



< Y^-k>t*-2 

k ~ 



Proof. Let t* > 2. In view of the second inequality in (J7.9|) 
(4.7) 



1 It* 
R < d n -i < — , so that q n < — < — . 

q n R R 



Next we compare the part of the graph of t i — >• ip r {t, v) that corresponds 
to tcosO > q n with the straight line T as defined on figure 5. They 
only differ when + kq n < tcos9 < q n -\ + (k + l)q n , and since ip r 
is a non increasing function, they differ by at most 



q n -i + kq n 
VA 7, — ,v 



((dn-i - (k - l)d n ) - R)q n . 



By (D, 



Hence 



cos 9 

< [d n _! - (k - l)d n ] - R<d n . 



t* 

< Ipr I — , V 

r 



Rq r , 



t* 

^ Qn ) "n 



< A 



cos# 



u 1 



<?n-l+fc9n<-7j <q„-i + (k + l)q n 



< Qnd-nXt* <(k+2)q n R < ^lfc>t*-2 > 



where the last inequality follows from ()4.Hj) and ()4.7|) . 

Next we estimate the difference 



I - R<i„ - [ p - <?n J 



1 - 



i? 



— t — 

d n -\ R 



Because of the second inequality in ()7.9|) . this difference is nonnegative. 
Because the map x t— > x + is a contraction, this difference is less than 



7? 



d 



Rq n + q n d n 



n-l 
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on the other hand both terms in the difference above vanish whenever 



R 



> 1n-i + (k + l)? n - Hence 



R 



< 



R 



< 



2 1 , 
+ - I It 



k+l^ k) "^<9n-l + (fc + lkn - k X 



k>t*-2 



where the penultimate inequality rests on (|4.3jl and the fact that R < 
d n -i. □ 



4.2. Step 2: applying the ergodic theorem. In this subsection 
again, we set r G (0, |) and pick G (0, |) such that a = tan# ^ 
Q; again we set v = (cos 9, sin 6 1 ) and i? = — ^ As in the previous 
subsection, we define n > 1 by the condition R E I n . In other words, 
we set n = N(a,R). With these assumptions, consider the expression 
Xr{%v) given by (@~5J), ie. 



R 



l-e A 



l(a,e *) _ ^* e -A (o,e 



with x = — In i? while Ao(a, e z ) and Ai(a, e x ) are defined as in (|3.7|) . 

Proposition 4.1. Lei iT G (0, 1). Lei f > 1; ffcera, /or ae. 9 G (0, f) 

swc/i that a = tan # ^ Q 



R" 



XRcosB T7; V 



hie /. /wlLUBl/ ViT 7 i2 



7T 



4- A 



1 / 1 + VI - z \ dz 

In . — V 1 — -2 -; 

1 - VT^z ) At* + z 

' z z \ dz 



1 - ^f\^~z J At* + z 



as e — > + , where v = (cos 9, sin 9). 
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Proof. The proof is based upon applying Proposition QQ] to the function 
/ defined by f(zi, z 2 ) = (1 — e z ' 2 — t*e~ Zl ) + , since, for each R* G (0, 1) 



1 f R * ft* \ dR 



lnel/. ^^U'^TT 



1 



lnel 



| lne| 

/(A (a,e- a! ),A 1 (a J e- a; ))dar. 



Starting from /, an elementary computation leads to F defined for each 
£ G (0, 1) as in Proposition 13. II bv 



F(0 = / /(I lnel - y, -y)dy = / (l - e"" - f e^+dy 

JO JO 

= £ / \c-i-nc*) + f- 



Assuming that t* > 1 and £ G (0,1), elementary computations show 
that 



C - 1 - t*£C 2 < for all C G R if 4t*f > 1 , otherwise 



,*t^2 ^ n :« 1 " VI - 4t*£ / , / 1 + VI - 4t*£ 



c _ i _ > o iff — ^— — ^ < c < 



2t*i ~ ~ 2t*i 



and hence 



i+V 1 - 4 **? 

2t*| /I I 



l-y/l-4i*g V t t 

2t*i v 



- - - 2 - t'Z | d£ 



1 - VI - 4t*£ 

+ 



1 + VI - 4t*£ 1 - VI - 4t*£ 
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12 
T 2 



12 

T 2 



4- ^ 



1 F(£)d£ 
1 + ^ 

-l/4t* 



Jo V 1 - VI - 4C« V 7 1+? 



l/4t* 



12 



4- 







i + vi - 4t*^ i - vi - 4t*£ y i + i 

dz 



2t*i 



1 / i n l±vO-7i 



1-vT 





4t* + z 



1 + VT^ 



□ 



4.3. Step 3: L 1 estimate of the remainder. The last ingredient in 
the proof of Theorem 11.11 consists in estimating the right hand side of 
(14. 6j) in average (integrating over the angle 9). 

Let R G (0, 1) and a G (0, 1) \ Q; we define n and k by the condition 
R £ I n ,h — below, this value of k is denoted by k(a, R). Equivalently, 

dN(a,R)( a ) ~ R 



n = N(pt,R) + l, k(a,R) 



dN(a,R)+l(&) 



+ 1. 



Lemma 4.3. Let m = m(6) G L°°((0, f )). Then, for each R G (0, 1) 
and each A > 1, one has 

-7T/4 



(4.9) 



lfc(tane,i?)>A m{6)d6 < 



2 m 



A 



Proof. The definition of k(a, R) implies in particular that 

dN(a,R)(tt) _ I 



fc(a,i2) < 



d 



N(a,R)+l 



a) 



1 



T N(a,R) a 



+ 1 



Hence, for each A > 1, one has 

dg — meas ({a G (0, 1) | k(a, R) > A}) 



< dg - meas fja G (0, 1) | < T N{a ' R) a < y~^j 



< dg — meas 



0. 



A- 1 



In 2 



In 



A 



A- 1 







as A — > +oo. Changing variables from a to 9 = arctana and using the 
classical inequality ln(l + z) < z leads to (|4.9jl . □ 
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4.4. Step 4: end of the proof. We conclude the proof of Theorem 
11.11 by bringing together the various ingredients described above. Let 

m e L°°([0, f]) such that m > and m(0)dO = 1. 

By Lemma 14.21 — and especially the inequality ()4.6|) there — and 
Lemma f4. 31 one has 

1 r e * r /4 ft* \ dr 

— J J i) r I -,(cos0,sin0)J m(9)d9— 

e * r /4 ft* \ dr 

Xr I — , (cos0, sin 9) I m(9)d9 



< 



In e J € J \r J r 

4 l re* rir/4 



< 



t* -2 | In 
Imlli 



1 f € r 1 dr 

J J 1 fc(tan 0, r / cos B)>t* -2?™(0) d9 — 



(t* - 3) 2 

By Proposition 14. II and dominated convergence, 



where 



1 — r/ / Xr (-,(cos0,sin0))m(0)d0 ► A(f 

lne| J t J \r J r 



1 - * J At* + z 

1 / z z \ dz 



as e — >• + . Hence 

i r* r /4 /r , .. . ..a dr . *ii//>ik- 



x r e n/^ ft* \ d T 

limsup- — r / / ip r — , (cos 0, sin 0) m{9)d9— < A(t*) + 



(t* - 3) 



while 



lim inf — ^ — p f r 7 % r f-,(cos0,sin0)^) m(9)d9— > A(t*)- ' S||, " I|l ' 



+o+ |lne|y e 7o \ r ' / r (t* — 3) 2 

As t* — > +oo, one has 

3 f l / 1 + v/T^ 
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Since the function t i— > ip r (t,v) is nonincreasing for all f G S* 1 and 
r G (0, |), one has, by using (|2.16|) . 

limsup — — r / / r I — , (cos 9, sin #) ) m(6)d8 — 

e^0+ I l ne l Je JO V r / r 

1 r 2s r /4 , A* „ \ dr 



limsup-- — r / / 4> r I — , (COS6 1 , sin^) m(9)d9- 
e _ + \\ne\J e Jo \r J 

1 /-2s /-Tr/4 ft* - L \ dr 

< limsup— — - / / ?/v I -, (cos 9, sin 9) I m[9)d9 — 

e ^ + |lne|y £ Jo \ r J r 

1 /- 2s W 4 ft* - s \ dr 

< limsup— — r / / ib r I , (cos 6 1 , sin #) m(#)<i6> — 

me J f In \ r r 



< A{t* -S) + 



\m\\ L o 



(t* - s - 3) 2 ' 

for each s G (0, |). Letting s — > + in the last inequality, one arrives 



at 



1 W4 / t * \ dr 

limsup— — r / / (j) r ( — , (cos 0, sin 0) I m{9)d9 



(4.10) 



lne| J(\ \ t J r 

<A(0 ' 8|H|i ~ 



(t* - 3) 2 ' 

A similar inequality holds for the liminf. By symmetry, the averaging 
in 9 in (J4.10)) can be done equivalently on (0,^) or in (0, 2tt), which 
eventually proves ()1.9|) . 

5. Applications to kinetic theory 

It has been proved in Theorem 2.1 of (see also [Hj) that the linear 
Boltzmann equation (ie. equation (10) of [T3]) does not govern the 
Boltzmann-Grad limit of the periodic Lorentz gas — unlike the case 
of a Lorentz gas with a random (Poisson) distribution of scatterers, 
where the linear Boltzmann equation was rigorously derived in jTT] 
and 5J. While Theorem 2.1 of is merely a negative result, we 
show below how to infer from Proposition 11.11 positive information on 
the asymptotic behavior of the periodic Lorentz gas in the Boltzmann- 
Grad limit. 

Define Q e = {ez \z G Z t }, and consider the transport equation 
d t f e + v ■ V x f e = , x G Q e , \v\ = 1 , 
(5.1) f e (t,x,v) = 0, x G <9fi £ , v-n x > 0, 

/ £ (0, X, V) = f n (x, V) iGfie, H=l. 
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Here, the unknown is f € = f e (t, x, v) while n x is the inward unit normal 
at point x G dQ t and f m is a given, nonnegative function of C C (R 2 x S 1 ). 
Physically, this is a variant of the periodic Lorentz gas where scatterers 
are replaced by holes (or traps) where impinging particles fall and thus 
are removed from the domain Q e . Obviously, for each t > 

(5-2) ||/e|k~,„ = ||/ in |k» • 

Reasoning as in [T3] suggests that /<= — > f in Lf x v weak-* where / 
solves the uniformly damped transport equation 

(5.3) dj + v VJ + f = on R; x R 2 x S 1 , f ]t=0 = f m , 

but this is ruled out by Theorem 2.1 of ^2]- Instead, Proposition 11.11 
suggests that the resulting damping rate should vanish in the limit as 
t — > +oo. This statement is made precise in the following theorem. 

Theorem 5.1. Let f in > belong to C C (R 2 x S 1 ) and let f e be, for 

each e 6 (0, V), the solution of \5.1)) . Then, for each nonnegative test 
function x £ Cl(R? x S 1 ), one has 

limsup // (- [ f r (t,x,v) — ] x( x 7 v)dxdv 

J J \ |lne| J e r J 

f(t, x, v)x(x, v)dxdv + O ( — 
limmf J J | j^-j [ f r (t,x,v dr 

f(t, x, v)x(x, v)dxdv + O ( — 





2f in (x-tv,v) 



as t — > +oo, where 

(5.4) f(t,x,v) 

n z t 

In particular, f satisfies 

(5.5) d t f + v-V x f + jf = 0, (t,x,v) E (0,+oo) x R 2 x S 1 
in the sense of distributions. 

Proof. First the solution of (|5.1j) is given by the formula 

(5.6) f e (t,x,v) = f m (x-tv,v)l Ti{x/e _ v) > t/t . 

In this formula, the exit time r e is considered as a function defined on 
Z e x S 1 with Z 2 -periodicity and extended by in x S 1 . 
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Let h = h(t, x, v) > belong to C°°(R + x R 2 x S 1 ), with support in 
R+ x [-L, L] 2 x S 1 ; 
(5.7) 

h(t,x, v)l Tr ^s. t - v )>*dxdv — / 4> r {-,v) / h(t,x,v)dxdv 



s 1 



< 



h(t, x, v)l 7 



+ 



/ez 2>/s JT 
7si V / ^ ez2 ■/ 



h(t, x, v)dx 



/ez 2 
4r»£ 



/ez 2 



\h(t,rl + ry,v) - h{t,rl,v)\l T , y _ v) >tdydv 

T 2 xS! 

| rZ + ry,t>) — h(t,rl,v)\(j) r ( ~,v ) oh/oh; 

<2r||V ;c /i|| LOO |5 1 |T 2 ^l [ _ LiL]2 (rO 

/ez 2 

<2r||V :E /i|| L oo| 1 S 1 |(L 2 + o(l)). 



By Theorem 11.11 



lim sup 



+0+ 



In el 



0r ( _ j v J / h(t,x,v)dxdv 



dr 



(5.1 



2 

8 



— // h(t,x,v)dxdv 



< 



t*-3 



with a similar estimate for the lim inf. 

Putting together (|5.7j) . (|5.8|) with the formula ()5.6j) establishes (|5.4jl . 

□ 

Theorem 15 . 1 1 can also be viewed as a result in homogenization. This 
remark leads to a comparison with the analogous situation — homog- 
enization of a diffusion process with Dirichlet boundary conditions in 
Z r — studied in (THj (in truth, the result obtained in [TH] is much 
more complete and satisfying than Theorem 15 .lp . Not surprisingly, the 
mathematical tools used in ^H] are of a very different nature than the 
ones in the present work. This however is by no means surprising and 
simply reflects the very different nature of the trajectories of a diffusion 
process and of those of a free transport equation. It is very likely that 
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the nature of the result in Theorem 15.11 — and its proof — would be 
deeply affected by adding some collision (ie. jump in velocity) process 
to the free transport between successive impingements on the obstacles. 

6. Final remarks and perspectives 

A first question left open in the present work is the existence of 
the limit ([1.7)1 . As the reader will have probably noticed in the proof 
of Theorem 11.1) the diameter r of the obstacle is treated as a time 
variable under renormalization — ie. transformation of the problem 
by the Gauss map. It may be that new insight concerning ([1.7)1 can be 
gained by using more specific properties of T than ergodicity as in the 
present paper. 

A second question is the existence of the limit in the sense of Cesaro, 
as commented upon in Remark 11.11 As can be seen from the proof of 
Theorem 11.11 proving this essentially amounts to being able to apply 
an ergodic theorem as in section [3] to functions of the form 

, ( d N (a, e) div(a,e)-i 

/ I £QN(a,e), " , " 

(In the present paper, ([4.3)1 essentially allows one to replace qu(a,e) 
by l/djv(a,e)-i)- This extension of the ergodic theorem in the present 
paper is postponed to a subsequent paper [5J and will be applied to 
the problem of the Lyapunov exponent for the Lorentz gas (see [5| for 
a presentation of this subject). 

Finally we would like comment on some very interesting, related 
work in [3| and jl]. 

The billiard problem considered in [2j is similar to studying the dis- 
tribution of free path lengths in Z r with the obstacle at the origin 
removed (ie. in Z r U 5(0, r/2)) for particles starting from the origin 
only. This is a quite different problem and thus the limit as r — > + of 
this distribution has nothing to do with the simple result in Theorem 
11.11 A related issue is studied in @]: for R > 0, consider the directions 
of lattice points (ie. point of Z 2 ) in the ball -6(0, R) of R 2 . These 
directions make a set of Nr angles 

< O < 0i,... < Nr = 2tt; 

the paper jlj computes the distribution of scaled differences of the 
form Nji(9k+i — Ok) as R — > +oo. Remarkably, these differences are 
far from being exponentially distributed — unlike in the case of angles 
picked at random in [0, 2ir). Consistently with the results in [S| and 
in [T2j, this observation stresses again the difference between the case 



24 



E. CAGLIOTI AND F. GOLSE 



of a random distribution of scatterers as studied in [TT] and that of a 
periodic distribution scatterers. 

7. Appendix: 
Background on continued fractions 

We recall below some basic facts and notations about continued frac- 
tions. Given a G (0, 1), the Gauss map is defined by 



(7.1) 



Ta = 

a 



it is known to be a uniquely ergodic transformation of (0, 1) with in- 
variant measure 

1 da 



(7.2) 



dg(a) 



In 2 1 + a 

The continued fraction expansion of a is 
(7.3) 

r ai, a 2 , a 3 , . . .] = with a k 



a 



ai + 



T k - l a 



k > 1 



a 2 + 



1 



a 3 + ... 

and the action of T is seen to correspond to the shift 

(7.4) a = [ai, a 2 , a 3 , . . .] h-> Ta = [a 2 , a 3 , a 4 , . . .] . 

The convergents of a are defined by the recursion formulas 

Qn+i = a n q n + Qn-i , % = , qi = 1 , 
p n+ i = a n p n + p n _i , p = 1 , pi = , 

and the corresponding error d n = \q n a — p n \ satisfies 

(7.6) d n+ i = -a n d n + d n -i , d = 1 , di = a . 
The first relation in ()7.5|) and ()7.6J) imply that 

(7.7) q n d n +i + q n +id n = 1 , n G N* . 

For each a G (0, 1) \ Q, the convergents of a are the sequence of best 
rational approximants of a. In other words 



(7.* 



\q n a -p n \= inf{ \qa - p\ | p, q G Z , < g < g n } . 



This implies in particular that the approximation by continued frac- 
tions is alternate in the sense that the algebraic value of the error is 

q n a -p n = (-ly^dn, n>0. 
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Also, for each a G (0, 1) \ Q, the sequence of errors satisfies the in- 
equalities 

(7.9) < d n < — . 

Qn + Qn+1 Qn+l 

The notations a n (a), p n {.ct), and d n (a) are used to emphasize the 

dependence of these quantities on a whenever necessary. 

Our discussion uses "renormalization" , ie. transforming a by the 
iterates of the Gauss map T. We have gathered some useful facts in 
the next lemma. 

Lemma 7.1. Let a G (0, 1); then 

• for each n G N, q n (Ta) = p n+ i(a) and d n+ i(a) = ad n {Ta); 

• for each nGN, 

d n (a) = \[T k a- 

k=0 

Proof. First observe that 

a n (Ta) = a n+ x(a) , n > 1 . 

This relation and (j7.5J) implies that the sequences q n (Ta) and p n+ i(a) 
satisfy the same recursion formulae; thus in order to check the first 
formula, it suffices to check it for both n = and n — 1 (in the latter 
case, one has qi(Ta) = [1/Ta] = a 2 = P2(®.)). The second formula is 
checked in the same way, and clearly implies the expression of d n (a) in 
terms of the T fc a's. □ 
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